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Abstract 

Molecular dynamics simulations are used to examine the relationship between water-like anoma- 
lies and the liquid-liquid critical point in a family of model fluids with multi-Gaussian, core-softened 
pair interactions. The core-softened pair interactions have two length scales, such that the longer 
length scale associated with a shallow, attractive well is kept constant while the shorter length scale 
associated with the repulsive shoulder is varied from an inflexion point to a minimum of progres- 
sively increasing depth. The maximum depth of the shoulder well is chosen so that the resulting 
potential reproduces the oxygen-oxygen radial distribution function of the ST4 model of water. As 
the shoulder well depth increases, the pressure required to form the high density liquid decreases 
and the temperature up to which the high-density liquid is stable increases, resulting in the shift 
of the liquid-liquid critical point to much lower pressures and higher temperatures. To understand 
the entropic effects associated with the changes in the interaction potential, the pair correlation 
entropy is computed to show that the excess entropy anomaly diminishes when the shoulder well 
depth increases. Excess entropy scaling of diffusivity in this class of fluids is demonstrated, showing 
that decreasing strength of the excess entropy anomaly with increasing shoulder depth results in the 
progressive loss of water-like thermodynamic, structural and transport anomalies. Instantaneous 
normal mode analysis was used to index the overall curvature distribution of the fluid and the 
fraction of imaginary frequency modes was shown to correlate well with the anomalous behaviour 
of the diffusivity and the pair correlation entropy. The results suggest in the case of core-softened 
potentials, in addition to the presence of two length scales, energetic and entropic effects associated 
with local minima and curvatures of the pair interaction play an important role in determining the 
presence of water-like anomalies and the liquid-liquid phase transition. 

PACS numbers: 64.70.Pf, 82.70.Dd, 83.10.Rs, 61.20.Ja 
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I. INTRODUCTION 



Water is characterized by well-known thermodynamic and kinetic liquid-state anomalies; 
for example, the rise in density on isobaric heating (density anomaly) and the increase in 
molecular mobility on isothermal compression (diffusivity anomaly) Since the anoma- 
lies of bulk water are connected with its behaviour as a solvent in chemical and biological 
systems, an understanding of the structural origins of such anomalous behaviour has at- 
tracted considerable attention^i^. While the anomalies of water were initially presumed to 
be uniquely connected to the hydrogen-bonded network of water-, there is now evidence 
that a number of hquids display water-like liquid-state anomalies, such as Te^ Ga, Bi^-, S^i^ 
GcisTegg^, silicajii"— , silicon^^ and BeF^^i'^^"'— . The generic relationships between struc- 
ture, entropy and mobility underlying this diverse set of liquids with water-like anomalies, 
can be understood in terms of the behaviour of the excess entropy (Sex), defined as the 
difference between the entropy [S) of the liquid and the corresponding ideal gas at the same 
density and temperature^^i^"— . A necessary condition for the fluid to show water-like ther- 
modynamic and transport anomalies is an excess entropy anomaly, corresponding to a rise 
in excess entropy on isothermal compression. The structural basis for the excess entropy 
anomaly is the existence of distinct forms of local order or length scales in the low- and 
high-density regimes; competition between the two types of local order results in a rise in 
excess entropy at intermediate densities. 

In addition to the singularity-free scenario for water-like thermodynamic and kinetic 
anomalies, it has been conjectured that the anomalies of water are due to the presence of 
a second liquid-liquid critical point, corresponding to the onset of a line of first-order phase 
transitions between high- and low-density phases of water—. The relationship between 
the liquid-liquid critical point and water-like anomalies can be addressed by considering 
minimal models of liquids with isotropic, pair-additive interactions that give rise to water- 
like anomalies, as well as liquid-liquid critical points^*^. Such models demonstrate that the 
presence of two repulsive length scales in the pair interaction is necessary in order to give rise 
to liquid-state anomalies. If, in addition to two repulsive length scales, the pair interaction 
has an attractive component, the fluid can show a liquid-liquid critical point (LLCP), in 
addition to the liquid-gas critical point. While two length scales in the pair interaction 
appears to be a necessary condition for seeing both the LLCP and water-like anomalies, it 
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is possible to design isotropic potentials with two length scales where appropriate variation 
of parameters can result in shifting either the LLCP or the water-like anomalies into the 
metastable or unstable regime^^i^i. 

In this paper, we study a family of liquids with continuous, core-softened pair interactions 
consisting of a hard core, a short-range shoulder and an attractive well at a larger separa- 
tion^. Since the potentials share a common functional form consisting of a sum of one 
Lennard- Jones and four Gaussian terms, we refer to them as the family of multi-Gaussian 
water-like liquids. By suitably varying the parameters, the outer attractive well can be left 
unchanged, while the shoulder can be progressively altered from being purely repulsive to a 
deep attractive well. As the shoulder shifts from being purely repulsive to more attractive, 
the anomalous regime in the pressure-temperature plane shrinks and disappears while the 
LLCP shifts to higher temperatures and lower pressures. The connection with atomistic 
models is made by ensuring that the limiting case of the double minimum potential that 
has no anomalies corresponds to an isotropic potential that reproduces the oxygen-oxygen 
radial distribution function of ST4 water^-. Using this set of anomalous fluids, we address 
a number of questions related to the features of the pair interaction, in addition to the two 
length scales, that control the temperature-pressure regime of the water-like liquid state 
anomalies versus the liquid-liquid critical point. Both the liquid-liquid critical point and the 
water-like anomalies require a change in the nature of local order in the liquid with density, 
and therefore two length scales in the case of core-softened fluids. The liquid-liquid critical 
point, however, depends on the energetic bias towards segregation of the two length scales 
with decreasing temperature. In contrast, the water-like liquid state anomalies require an 
excess entropy anomaly, involving a continuous transformation of the liquid from low- to 
high-density through a range of quasi-binary states reflecting the competition between two 
length scales in the intermediate regime^. In order to understand the relationship between 
the interaction potential, the water-like liquid state anomalies and the liquid-liquid critical 
point, it is therefore necessary to consider the temperature-dependent stabilization of the 
low- and high-density length scales as well as the density-dependent changes in the entropy 
of the system. 

In order to understand how the change in interaction potential within the multi-Gaussian 
water models affects the thermodynamic and kinetic water-like anomalies, it is necessary 
to map out the excess entropy anomaly for the different model fluids. We use the pair 
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correlation entropy as a simple structural estimator of the excess entropy, defining it for a 
one-component fluid of structureless particles as: 

= -27rp J {g{r) \ng{r) - [g{r) - l]}r'^dr (1) 

where g{r) is the radial distribution function. It is typically the dominant contribution to 
the excess entropy of a fluid expressed as a multi-particle correlation expansion of the form: 

Sex = S — Sid = S2 + S3 + . . . (2) 

where Sn is the entropy contribution due to n-particle spatial correlations^^"— . The excess 
entropy and mobility anomalies are linked by excess entropy scaling relations of the form 

X* = Aexp{a{S,jNkB)) , (3) 

where X* are dimensionless transport properties with either macroscopic (Rosenfeld) or 
microscopic (Dzugutov) reduction parameters and the scaling parameters, a and A, depend 
on the functional form of the underlying interactions22."ii. 

As an additional means to relate the interaction potential to the liquid-state properties, 
we characterize the potential energy surface (PES) of the multi-Gaussian family of water- 
like liquids using instantaneous normal mode analysis. In the Instantaneous Normal Mode 
(INM) approach, the key quantity is the ensemble-averaged curvature distribution of the 
PES sampled by the system. For a system of particles, the mass-weighted Hessian as- 
sociated with each instantaneous configuration is diagonalized to yield 3A^ normal mode 
eigenvalues and eigenvectors and the ensemble-average of this distribution is referred to as 
the INM spectrum. The INM spectrum of a liquid will have a substantial fraction of unsta- 
ble modes with negative eigenvalues that simulations suggest is strongly correlated with the 
diffusivity^"— . Random energy models of liquids also suggest that for supercooled liquids 
there will be a connection between the fraction of imaginary modes, the diffusivity and the 
configurational entropy^"^. Our previous work on INM analysis of a core-softened water-like 
fluid demonstrated that the instantaneous normal mode spectra carry significant informa- 
tion on the dynamical consequences of the interplay between length scales characteristic of 
anomalous fluids^^. We have therefore performed an INM analysis of the multi-Gaussian 
water-like fluids to understand the relationship between the interaction potential, anomalies 
and the liquid-liquid critical point. 



The paper is organized as follows. The computational details of the simulations and the 
equation of state data for the multi-Gaussian family of water models are summarized in 
Section II. Section III contains the results and the conclusions in Section IV. 



II. THE MODEL 



A. The potential 



The multi- Gaussian family of water-like fluids is defined by pair-additive, continuous, 
core-softened interactions with the functional form 



f/(r) 



hj exp 



r — Cn 



(4) 



8 



* 

U 4 


1.2 1.8 2.4 3.0 
r 

FIG. 1: Interaction potential obtained by changing parameters hi in Eq. @. The potential and 
the distances are in dimensionless units U* = U/'j and r* = r/rg. 

The first term is a Lennard- Jones potential-like and the second one is composed of four 
Gaussians, each of width wj centered at Cj. The potential and the distances are given in 
dimensionless units, U* = U/'y and r* = t/tq where 7 is the energy scale and tq is the 
length scale chosen so the closest approach between particles is about r* = 1, i.e., so that 
the second length scale associated with the repulsive shoulder remains the same. Here we 
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TABLE I: Parameters hi for potentials A, B, C and D. 



Potential 


Value oi hi 


A 


0.25 h\'f 


B 


0.50 /i^^ 


C 


Q.15h\'f 


D 


1.00 /i^^ 



TABLE II: Parameters for potentials A, B, C and D in units of A and of kcal/mol. 



Parameter 


Value 


Parameter 


Value 


a 


9.056 


Wl 


0.253 


b 


4.044 


W2 


1.767 


e 


0.006 


W3 


2.363 


a 


4.218 




0.614 


Cl 


2.849 


hl'f 


-1.137 


C2 


1.514 


h2 


3.626 


C3 


4.569 


hs 


-0.451 


C4 


5.518 


/14 


0.230 



use e/7 = 0.02 and cr/ro = 1.47. Modifying hi in the Eq. (jlj) allows us to change the depth 
of the hard-core well, as illustrated in Fig. [1] while keeping the shape and location of the 
attractive well constant. We report here results for four different values for hi and they are 
expressed as a multiple of a reference value h]^^ as shown in the Table [H For all the four 
cases the values of a, b, {cj, Wj} with j = 1, . . . , 4 and /i^''^'^. Table [III gives the parameter 
values in A and kcal/mol consistent with reproducing the oxygen-oxygen radial distribution 
of ST4 water using case 
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B. The simulation details 



The properties of the system were obtained by NVT molecular dynamics using Nose- 
Hoover heat-bath with coupling parameter Q — 2. The system is characterized by 500 
particles in a cubic box with periodic boundary conditions, interacting with the intermolec- 
ular potential described above. All physical quantities are expressed in reduced units and 
defined as 

i(m/7)V2 



P 



To 

7 
pro 



7 

* 3 

P = prt 

D' = ^ . (5) 

in 

Standard periodic boundary conditions together with predictor-corrector algorithm were 
used to integrate the equations of motion with a time step At* = 0.002 and potential cut 
ofi' radius r* = 3.5. The initial configuration is set on solid or liquid state and, in both 
cases, the equilibrium state was reached after t*^ = 1000 (what is in fact 500000 steps since 
At* — 0.002) . From this time on the physical quantities were stored in intervals of Ai^ — 1 
during t*j^ — 1000. The system is uncorrelated after t*^ — 10, as judged from the velocity 
auto-correlation function. 50 descorrelated samples were used to get the average of the 
physical quantities. 

At each state point, 100 configurations were sampled and used to construct the instanta- 
neous normal mode spectra and associated quantities. We repeated the calculation for some 
state points using 500 configurations and found no significant difference. 

C. Instantaneous Normal Modes Analysis 

The potential energy of configuration r near Fq can be written as a Taylor expansion of 
the form: 

U{r) = C/(ro) - F • z + ^r"^ • H • z (6) 
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where Zj = y'm7(ri — ro) are the mass-scaled position coordinates of a particle i. The 
first and second derivatives of U{r) with respect to the vector z are the force and the 
Hessian matrix, denoted by F and H respectively. The eigenvalues of the Hessian H are 
{{uf},i = 1,3A^) representing the squares of normal mode frequencies, and W(r) are the 
corresponding eigenvectors. In a stable solid, fq can be conveniently taken as the global 
minimum of the potential energy surface U{R), which implies that F = and H has only 
positive eigenvalues corresponding to oscillatory modes. The INM approach for liquids 
interprets r as the configuration at time t relative to the configuration fq at time to- Since 
typical configurations, tq are extremely unlikely to be local minima, therefore F 7^ and 
H will have negative eigenvalues. The negative eigenvalue modes are those which sample 
negative curvature regions of the PES, including barrier crossing modes. The ensemble- 
averaged INM spectrum, (/(a;)), is defined as 

Quantities that are convenient for characterizing the instantaneous normal mode spectrum 
are: (i) the fraction of imaginary frequencies (Fi), defined as 



F,= / f{io)duj (8) 

J im 

and the Einstein frequency (uje), given by 

= J uj'^f(uj)duj 

- <^^"> (9) 



m{3N - 3) 

where the integral is performed over the entire range of frequencies, real as well as imaginary. 



III. RESULTS 



A. Phase Diagram 

Figure [2] illustrates the pressure-temperature phase diagram for the four cases of the 
potentiaF^. Due to the presence of the attractive interaction, all four cases have a liquid- 
gas transition with an associated critical point that is not shown here. In addition, all the 
four model liquids studied here have a liquid-liquid critical point. Cases A, B and C have 
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water-like density and diffusional anomalies. The solid bold lines represent the locus of 
temperatures of maximum density (TMD) for different isobars. State points enclosed by 
the TMD locus represent the regime of density anomaly within which {dp/dT)p > 0. The 
maximum temperature along the TMD locus, denoted by T^^fj-,, is the threshold temperature 
for onset of the density anomaly. The dot-dashed lines are the temperatures of maximum 
and minimum diffusivity along different isotherms^. 

This overall change in the nature of the liquid-state phase diagram for the four multi- 
Gaussian liquids is summarized in Figure |3l Clearly, as the second length scale shifts from 
an inflexion point on the repulsive shoulder to a well with progressively increasing depth 
and curvature, the region of liquid state anomalies shrinks and disappears. The figure 
illustrates how the pressure and temperature associated the liquid-gas and liquid-liquid 
critical points vary with the potentials A,B,C and D. The same graph also shows that as 
the shoulder becomes deeper, the maximum temperature of the TMD locus, which marks 
the onset temperature for thermodynamically anomalous behaviour, approaches the liquid- 
liquid critical point. 

Since the thermodynamic and mobility anomalies of water are correlated, we first focus 
on understanding the thermodynamic condition for the presence of density anomaly. This 
may be stated as: 

^ = -^>0 (10) 

dp NKt ^ ' 

where a is the thermal expansion coefficient and Kt is the isothermal compressibility. For 
the system to have a large anomalous region, the ratio u/Kt should be therefore large and 
negative. Near the critical point, the compressibility, Kt-, and thermal expansion coefficient, 
aTi diverge, however the compressibility diverges with a large exponent making the ratio 
zero. In this case, the condition given by Eq. (fTOj) can not be fulfilled. This suggests that 
near the liquid-liquid critical point the system prefers to undergo a phase separation into 
high- and low-density liquids, rather than show a smooth entropy anomaly. 



B. Excess Entropy and Pair Correlations 

As discussed in the previous section, the density anomaly corresponds to a set of state 
points for which {dS/dp)T > 0. The total entropy is a sum of the ideal (Sid) and excess 
(Sex) contributions. Since Sid decreases monotonically with increasing density, therefore a 
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FIG. 2: Pressure-temperature phase diagram for cases A, B, C and D. The thin sohd hues are the 
isochores 0.30 < p* < 0.65. The liquid-hquid critical point is the dot, the locus of temperatures of 
maximum density is the solid thick line and the locus of diffusion extrema is the dot-dashed line. 



density anomaly must imply the presence of an excess entropy anomaly, {dSex/ dp)T > 0. 
Errington et al. have further shown that the strength of the excess entropy anomaly required 
to give rise to density anomaly is given by the condition = id{Scx/NkB)/d\np)j, > 1— . 
By approximating the excess entropy with the two-body correlation contribution of S2 [see 
Eq. (1)], we relate the structural information in the radial distribution function of the fluid 
to the thermodynamic behaviour. 

Figure m illustrates the S2(p) = S2/NkB versus p* for various temperatures and for the 
potentials A, B, C and D. For cases A, B and C, at low temperatures, there is a rise in 
excess entropy on isothermal compression characteristic of water-like liquids^^'i^i^i^i^ that 
contrasts with the behaviour of simple liquids where free volume arguments are sufficient to 
justify a monotonic decrease in entropy on isothermal compression. For case D, no anomaly 



11 




0.5 1.0 1.5 2.0 2.5 

* 

T 

FIG. 3: Pressure versus temperature locations of the critical points for the potentials A — D. 

is observed in the pair entropy even at low temperatures. The progressive attenuation of 
the anomalies on going from case A to case D, is illustrated in Figure |5] which compares the 
behavior of the pair entropy versus density for all studied potentials at a given temperature, 
T* = 0.9. This graph together with Figure [3] indicates that as the maximum temperature 
at the TMD line approaches the liquid-liquid critical temperature, the pair entropy curve 
becomes more flat and the anomalous behavior disappears. 

The origin of the pair entropy anomaly in fluids with two length scales can be explained 
in terms of a competition between two length scales at intermediate densities. Only a single 
length scale dominates in the low- and high-density limits while at intermediate densities, 
where both length scales are present, can be regarded as quasi-binary systems with a mixing 
entropy. The radial distribution functions shown in our previous study clearly demonstrate 
the presence of two length scales. They also show that with increasing temperature, the 
shorter length scale peak of g{r) becomes more prominent in cases A, B and C. In contrast, 
in case D, both length scales associated with the first and second peak of the g{r) broaden 
with increasing temperature as a consequence of which there is no emergence of an anomaly 
with decreasing temperature. 

The crucial question to ask in the multi-Gaussian family of water models is why, de- 
spite the presence of two length scales at intermediate densities, the pair entropy anomaly 
is progressively lost as the shoulder goes from being an inflexion point to a minimum with 
about twice the depth as the outer, attractive well. Clearly the rise in entropy with isother- 
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mal compression due to mixing of two length scales is counteracted by additional effects. 
To understand this we note that the entropy of an one- dimensional harmonic oscillator of 
frequency uj is given by: 



The increasing curvature of the short-range minimum, relative to the attractive minimum, 
implies that a pair separated trapped in the shoulder minimum will have lower vibrational 
entropy than one trapped in the broad shallow attractive minimum. As a consequence, at 
intermediate densities, while the presence of two length scales will increase entropy, the loss of 
entropy when the pairs are located in the short-range minimum will tend to decrease entropy. 
As the shoulder minimum becomes deeper, the second effect becomes more important and 
the excess entropy anomaly disappears. In systems such as the two-scale linear ramp, such 
curvature-dependent effects will be absent. 

It is also interesting to consider the shifting of the liquid-liquid critical point to lower 
pressures and higher temperatures. For a temperature-driven phase separation into low- 
density liquid (LDL) and high-density liquid (HDL), increasing energetic stabilization of 
one length scale relative to the other is required. In case A, this is clearly due to the outer 
attractive well with a depth of about At/* 0.3 and T* ^ 0.3. In case D, this is due to 
the shorter length scale with a well depth of about A* = 1.00 and T* ^ 0.8. For a shallow 
shoulder, the high density liquid is stabilized under high pressure. Within the HDL phase 
particles are occupying the shoulder scale. The density anomalous region, characterized by 
having particles in the two scales occurs at the pressure range of the low density liquid phase. 
As the shoulder scale becomes deeper, less pressure is needed to form the high density liquid. 
The LDL phase occupies a smaller pressure range and therefore the density anomalous region 
shrinks. For a very deep shoulder well as in case D, the HDL requires no pressure to be 
formed, the LDL is at negative pressures and the anomalous density regime disappears. 

C. Diffusion and Rosenfeld Reduction Parameter 

Previously we have shown that the diffusion coefficient in the cases A, B and C decreases 
with the decrease of the density for a certain range of densities-^ . The region in the pressure 
temperature phase diagram limited by the maxima and minima of the diffusion coefficient 
is illustrated as dot-dashed lines in Fig. O In the case D the diffusion coefficient increases 




(11) 
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FIG. 4: Pair entropy versus density for for the cases A, B, C and D for various temperatures. 

with the decrease of the density as in normal hquids. It is interesting to notice that the same 
behavior is also observed in the pair entropy suggesting that the anomalies present in these 
two quantities might be related. In order to check this hypothesis we now consider the scaling 
relationship between the diffusivity and the pair entropy. Using the Rosenfeld macroscopic 
reduction parameters for the length as p"^/^ and the thermal velocity as (fc^T/m)^/^, the 
dimensionless diffusivity is defined as 

The scaling of the reduced diffusivity, D^i with pair entropy, is illustrated in Fig. Ei 
Previous results for core-softened fluids^ suggest that /S.S = Sex — S2 tends to be density 
dependent in anomalous fluids, resulting in a stronger isochores dependence when \ia{Dfj) is 
plotted against 5*2, rather than against Sex- In the present study, we have computed only 5*2 
and therefore Fig. [6] shows scaling with respect to 5*2. For case A, the collapse of data from 
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FIG. 5: Pair entropy versus density for the cases A, B, C and D at T* = 0.9. 

all the state points on a single line is quite good. As we progress from case A to case D, the 
isochores dependence of the scaling parameters becomes more pronounced suggesting that 
the density dependence of AS* increases on going from case A to case D. 

D. The Instantaneous Normal Mode Spectrum 

The variation in anomalous behaviour in the multi-Gaussian family of water-like liquids 
studied here suggests that in addition to length scales, we need to look at other features of 
the pair potential e.g. its first and second derivatives. Instantaneous normal mode analysis 
provides a way to summarize information on the curvature distribution of the potential 
energy landscape. In Figure [TJ we show the INM spectra of liquids bound by the four 
potentials (A, B, C and D) at a common state point of p* = 0.50 and temperature T* = 0.8. 
The crucial features are as follows: 

• A low-frequency split peak in the real branch centered at about u = 10, that does not 
vary significantly between the four cases and must reflect modes associated with the 
outer attractive well; 
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FIG. 6: Diffusion in Rosenfeld units as a function of — So for the cases A, B, C and D. 



• A high-frequency peak in the real branch, centered at approximately 30, 35, 40 and 
50 for cases A, B, C and D respectively, which must correspond to motion in neigh- 
borhood of the shoulder length scale. As the curvature of the short-range minimum 
increases, this features shifts to higher frequencies and becomes more prominent; 

• The imaginary branch reflects regions of negative curvature in the neighborhood of 
barriers and inflexion points. Case A, where there is no barrier in the pair interaction 
but only an inflexion point has a single peak like a simple liquid. However, this peak 
is broad because of the core-softened repulsive wall and the fraction of imaginary 
modes is large. As the barrier between the short and long length scales becomes more 
pronounced in the pair interaction, the second peak in the imaginary branch becomes 
more prominent. 

Thus the real branch is dominated by vibrational modes associated with motion in the 
attractive and shoulder length scales while the imaginary modes branch is dominated by 
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negative curvature modes associated with transitions between the shoulder and attractive 
wells. For the case A, the real branch has three peaks related to the three basins: the 
shoulder scale, the attractive scale and a second attractive scale located at further distance 
in Fig. [TJ It has just one imaginary peak that indicating that transitions between the two 
length scales do not require local barrier crossing. For the cases B and C the imaginary 
branch has two peaks suggestive of modes connecting between the shoulder scale, attractive 
scale and second attractive scales. The peak with largest frequency in the real branch has 
larger frequency in the case C than in the case A and is related to the shoulder scale. For 
the case D, the shoulder is deep and so the frequency related to the shoulder scale has a 
very high frequency. The imaginary branch has two distinct oscillation modes that exclude 
transitions between the shoulder scale and the other scales and therefore no anomalies are 
expected. 

The above discussion suggest that INM spectra carry fairly detailed information on the 
dynamics of transitions between the two length scales. The two features which are a com- 
pact signature of INM spectra are the Einstein frequency and the fraction of imaginary 
modes. Isotherms of the Einstein frequency as a function of density for all four cases show 
a monotonic increase with density and do not show any significant signatures of the water- 
like anomalies. The fraction of imaginary modes, in contrast, correlates strongly with the 
anomalous behaviour of the pair entropy and the diffusivity. Figure [8] shows the Fi curves 
versus density for various isotherms of all four multi-Gaussian model fluids studied here. The 
parallel behaviour of the S2{p) and Fi{p) curves at corresponding isotherms is immediately 
obvious, though the Fi{p) have a stronger non- monotonic behaviour than S2{p) curves. This 
can be seen most clearly for a high-temperature isotherm. 

IV. CONCLUSIONS 

This paper examines the relationship between water-like anomalies and the liquid-liquid 
critical point in a family of model fluids with multi-Gaussian, core-softened pair interactions. 
The pair interaction in this family of liquids is composed of a sum of Lennard- Jones and 
Gaussian terms, in such a manner that the longer length scale associated with a shallow, 
attractive well is kept constant while the shorter length scale associated with the repulsive 
shoulder changes from an inflexion point to a minimum of progressively increasing depth. 
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FIG. 7: Normal models versus frequency for the four studies cases. The density is fixed, p* = 0.50 
in all the cases and the temperature is varied. 



The maximum depth of the shoulder length scale is chosen so that the resulting potential 
reproduces the oxygen-oxygen radial distribution function of the ST4 model of water. As the 
energetic stabilization of the shoulder length scale increases, the liquid-liquid critical point 
shifts to higher temperatures and lower pressures. Simultaneously, the temperature for onset 
of the density anomaly decreases and the region of liquid state anomalies in the pressure- 
temperature plane diminishes. The condition for the presence of anomalies is inconsistent 
with divergences near a critical point, so that in the limiting case of maximum shoulder well 
depth, the anomalies disappear. 

To understand our results for the phase diagram and liquid-state anomalies of the multi- 
Gaussian family of water-like fluids, it is important to note that, in addition to the presence 



18 



0.25 
0.20 
F.0.15 
0.10 
0.05 

" 0.3 0.4 0.5 0.6 " 0.3 0.4 0.5 0.6 

P P 

0.25 
0.20 
FjO.lS 
0.10 
0.05 

^ 0.3 0.4 0.5 0.6 " 0.3 0.4 0.5 0.6 

P P 

FIG. 8: Fraction of imaginary modes versus density for fixed temperatures, T* = 
0.20,0.30,0.40,0.50,0.60,0.70,0.80,0.90 and 1.10 from bottom to top for cases A,B and C. For 
case D the start at T* = 0.50. 

of two length scales, it is necessary to consider the energetic and entropic effects as de- 
termined by local minima and curvatures of the pair interaction. As the shoulder depth 
increases, the pressure required to form the high density liquid decreases and the temper- 
ature up to which the high-density liquid is stable increases. This explains the shift of the 
liquid-liquid critical point to much lower pressures and higher temperatures. To understand 
the entropic effects associated with the changes in the interaction potential, we computed the 
pair correlation entropy and demonstrated the attenuation of the excess entropy anomaly as 
the shoulder length scale changed from an inflexion point to a deep minimum. In conjunc- 
tion with Rosenfeld-scaling of transport properties, this is consistent with the progressive 
loss of water-like thermodynamic, structural and transport anomalies. The excess entropy 
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anomaly in two-scale, isotropic fluids is due to a rise in entropy as a result of competition 
between two length scales at intermediate densities. In the case of continuous potentials, the 
vibrational entropy associated with the two length scales becomes important. To index the 
overall curvature distribution in the liquid, we have used instantaneous normal mode analy- 
sis and shown the fraction of imaginary frequency modes correlates well with the anomalous 
behaviour of the diffusivity and the pair correlation entropy. A detailed analysis of the 
INM spectrum shows that as the shoulder well increases in depth, there is a simultaneous 
rise in the positive curvature associated with the shoulder minimum as well as the negative 
curvature of the barrier separating the shoulder minimum from the attractive minimum. 
Consequently, the vibrational entropy associated with pairs of particles separated by the 
shoulder distance decreases, relative to that of pairs trapped in the outer attractive well. 
Therefore the mixing entropy due to the presence of two length scales is counteracted by 
the changes in vibrational entropy associated by the two length scales. 

A general conclusion that emerges from this study is that even though the ratio between 
the two length scales is important for locating the temperature range of the anomalies^, 
additional energetic and entropic effects associated with local minima and curvatures of the 
pair interaction can play an important role. The liquid-liquid phase separation depends on 
the relative energies associated with the two length scales whereas the water-like anomalies 
depend upon a continuous rise in entropy as a function of isothermal compression. A number 
of recent studies of core-softened fluids illustrate this conclusion. For example, energetic 
and entropic effects play a very different role in the discrete and discontinuous versions 
of the shouldered well potential^^. In the discrete case, the enthalpic implications do not 
change significantly and the liquid-liquid critical point is not significantly different in the two 
systems. In contrast, the continuous potential allows for a smooth transformation through 
a range of quasi-binary states from low- to high-density and shows water-like anomalies. A 
more recent study of core-softened fiuids^^ shows that increasing the depth of the attractive 
well, while leaving the shoulder feature constant, results in disappearance of the anomalies 
while shifting the liquid-liquid critical point to lower pressures and higher temperatures. 
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